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A novel rf-SQUID flux qubit that is robust against fabrication variations in Josephson junction 
critical currents and device inductance has been implemented. Measurements of the persistent 
current and of the tunneling energy between the two lowest lying states, both in the coherent and 
incoherent regime, are presented. These experimental results are shown to be in agreement with 
predictions of a quantum mechanical Hamiltonian whose parameters were independently calibrated, 
thus justifying the identification of this device as a flux qubit. In addition, measurements of the 
flux and critical current noise spectral densities arc presented that indicate that these devices with 
Nb wiring are comparable to the best Al wiring rf-SQUIDs reported in the literature thusfar, with 
a 1// flux noise spectral density at 1 Hz of I.SIq's Z^'&o/v'Hz. An explicit formula for converting the 
observed flux noise spectral density into a free induction decay time for a flux qubit biased to its 
optimal point and operated in the energy eigenbasis is presented. 



I. MOTIVATION 

Experimental efforts to develop useful solid state quan- 
tum information processors have encountered a host 
of practical problems that have substantially limited 
progress. While the desire to reduce noise in solid state 
qubits appears to be the key factor that drives much of 
the recent work in this field, it must be acknowledged that 
there arc formidable challenges related to architecture, 
circuit density, fabrication variation, calibration and con- 
trol that also deserve attention. For example, a qubit 
that is inherently exponentially sensitive to fabrication 
variations with no recourse for in-situ correction holds 
little promise in any large scale architecture, even with 
the best of modern fabrication facilities. Thus, a qubit 
designed in the absence of information concerning its ul- 
timate use in a larger scale system may prove to be of 
little utility in the future. In what follows, we present an 
experimental demonstration of a novel superconducting 
flux qubit ^ that has been specifically designed to address 
several issues that pertain to the implementation of a 
large scale quantum information processor. While noise 
is not the central focus of this article, we nonetheless 
present experimental evidence that, despite its physical 
size and relative complexity, the observed flux noise in 
this flux qubit is comparable to the quietest such devices 
reported upon in the literature to date. 

It has been well established that rf-SQUIDs can be 
used as qubits given an appropriate choice of device pa- 
rameters. Such devices can be operated as a flux biased 
phase qubit using two intrawell energy levels^ or as a 
flux qubit using any pair of interwell levels^. This ar- 
ticle will focus upon an experimental demonstration of 
a novel rf-SQUID flux qubit that can be tuned in-situ 
using solely static flitx biases to compensate for fabrica- 
tion variations in device parameters, both within single 
qubits and between multiple qubits. It is stressed that 



this latter issue is of critical importance in the develop- 
ment of useful large scale quantum information proces- 
sors that could foreseeably involve thousands of qubits^. 
Note that in this regard, the ion trap approach to build- 
ing a quantum information processor has a considerable 
advantage in that the qubits are intrinsically identical, 
albeit the challenge is then to characterize and control 
the trapping potential with high fidelity*. While our re- 
search group's express interest is in the development of a 
large scale superconducting adiabatic quantum optimiza- 
tion [AQO] processor^'^, it should be noted that many 
of the practical problems confronted herein are also of 
concern to those interested in implementing gate model 
quantum computation [GMQC] processors'" using super- 
conducting technologies. 

This article is organized as follows: In Section II, a 
theoretical argument is presented to justify the rf-SQUID 
design that has been implemented. It is shown that this 
design is robust against fabrication variations in Joseph- 
son junction critical current. Second, it is argued why it 
is necessary to include a tunable inductance in the flux 
qubit to account for differences in inductance between 
qubits in a multi-qubit architecture and to compensate 
for changes in qubit inductance during operation. There- 
after, the focus of the article shifts towards an experimen- 
tal demonstration of the rf-SQUID flux qubit. The archi- 
tecture of the experimental device and its operation are 
discussed in Section III and then a series of experiments 
to characterize the rf-SQUID and to highlight its control 
are presented in Section IV. Section V contains measure- 
ments of properties that indicate that this more complex 
rf-SQUID is indeed a flux qubit. Flux and critical cur- 
rent noise measurements and a formula for converting the 
measured flux noise spectral density into a free induction 
(Ramsey) decay time are presented in Section VI. A sum- 
mary of key conclusions is provided in Section VII. De- 
tailed calculations of rf-SQUID Hamiltonians have been 
placed in the appendices. 
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II. RF-SQUID FLUX QUBIT DESIGN 

The behavior of most superconducting devices is gov- 
erned by three types of macroscopic parameters: the 
critical currents of any Josephson junctions, the net ca- 
pacitance across the junctions and the inductance of the 
superconducting wiring. The Hamiltonian for many of 
these devices can gcncrically be written as 
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where Ci, Eji = /i(f>o/27r and Ii denote the capaci- 
tance, Josephson energy and critical current of Joseph- 
son junction i, respectively. The terms in the first sum 
are readily recognized as being the Hamiltonians of the 
individual junctions for which the quantum mechani- 
cal phase across the junction ipi and the charge col- 
lected on the junction Qi obey the commutation relation 
[<i>o</5i/27r, Qj] = ihSij. The index n in the second sum- 
mation is over closed inductive loops. External fluxes 
threading each closed loop, <i>Jj, have been represented as 
phases ip^-^ = 27r<i>^/<i>o- The quantum mechanical phase 
drop experienced by the superconducting condensate cir- 
culating around any closed loop is denoted as ipn- The 
overall potential energy scale factor for each closed loop 
is given by t/„ = (<i>o/27r)^/L„. Here, L„ can be either 
a geometric inductance from wiring or Josephson induc- 
tance from large junctions*. Hamiltonian (1) will be used 
as the progenitor for all device Hamiltonians that follow. 



Compound-Compound Josephson Junction 
Structure 



A sequence of rf-SQUID architectures are depicted in 
Fig. 1. The most primitive version of such a device is 
depicted in Fig. la, and more complex variants in Figs.l 
b and Ic. For the single junction rf-SQUID (Fig. la), the 
phase across the junction can be equated to the phase 
drop across the body of the rf-SQUID: ipi = ipq. The 
Hamiltonian for this device can then be written as 
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with the qubit inductance Lg = Lbodyi qubit capacitance 
Cq = Ci and qubit critical current 1^ = Ii in this particu- 
lar case. If this device has been designed such that /3 > 1 
and is flux biased such that t/?^ « tt, then the poten- 
tial energy V{(fiq) will be bistable. With increasing f3 an 
appreciable potential energy barrier forms between the 
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FIG. 1: (color online) a) A single junction rf-SQUID qubit. 
b) Compound Josephson Junction (CJJ) rf-SQUID qubit. c) 
Compound-Compound Josephson Junction (CCJJ) rf-SQUID 
qubit. Junction critical currents U and junction phases Lpi 
(1 < i < 4) as noted. Net device phases are denoted as ifa, 
where a £ {£,r,q). External fluxes, cO^^ are represented as 
phases tp^ = 27r"I>^/"l>o, where n G (L, _R, cjj, ccjj, q). Induc- 
tance of the rf-SQUID body, CJJ loop and CCJJ loop are 
denoted as Lbody, icjj and I/ccjj, respectively. 



two local minima of V{ipq), through which the two low- 
est lying states of the rf-SQUID may couple via quantum 
tunneling. It is these two lowest lying states, which are 
separated from all other rf-SQUID states by an energy of 
order of the rf-SQUID plasma energy hujp = h/ yjLqC\, 
that form the basis of a qubit. One can write an effective 
low energy version of Hamiltonian (2a) as^ 



^9 = ^2 t^^^ 



(3) 



where e = 2 (<i>;^ - <i>o/2), is the magnitude of 
the persistent current that flows about the inductive q 
loop when the device is biased hard [e » A,] to one side 
and Ag represents the tunneling energy between the oth- 
erwise degenerate counter-circulating persistent current 
states at <i>|^ = <i>o/2. 

A depiction of the one-dimensional potential energy 
and the two lowest energy states of an rf-SQUID at de- 
generacy (<i>|^ = <i>o/2) for nominal device parameters is 
shown in Fig. 2. In this diagram, the ground and first ex- 
cited state are denoted by \g) and |e), respectively. These 
two energy levels constitute the energy eigenbasis of a 
fiux qubit. An alternate representation of these states, 
which is frequently referred to as either the fiux or persis- 
tent current basis, can be formed by taking the symmet- 
ric and antisymmetric combinations of the energy eigen- 
states: |i) = (|g) + |e)) /V2 and |T) = (|<?) - |e)) / A 
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FIG. 2: (color online) Depiction of the two lowest lying states 
of an rf-SQUID at degeneracy (e = 0) with nomenclature for 
the energy basis (|<7),|e)) and flux basis (|i),|T)) as indicated. 



which yield two roughly gaussian shaped wavefunctions 
that are centered about each of the wells shown in Fig. 2. 
The magnitude of the persistent current used in Eq. (3) 
is then defined by = |(T| - $o/2) /L, |T)|. The 



tunneling energy is given by 



{9\-n,\9). 



The aforementioned dual representation of the states 
of a flux qubit allows two distinct modes of operation of 
the flux qubit as a binary logical element with a logical 
basis deflned by the states |0) and In the first mode, 
the logical basis is mapped onto the energy eigenbasis: 
|0) \g) and |1) — + |e). This mode is useful for optimiz- 
ing the coherence times of flux qubits as the dispersion 
of Hamiltonian (3) is flat as a function of <I>g to flrst or- 
der for e w 0, thus providing some protection from the 
effects of low frequency flux noise^". However, this is not 
a convenient mode of operation for implementing inter- 
actions between flux qubits^^'^^. In the second mode, 
the logical basis is mapped onto the persistent current 
basis: |0) \[) and |1) This mode of opera- 

tion facilitates the implementation of inter-qubit interac- 
tions via inductive couplings, but does so at the expense 
of coherence times. GMQC schemes exist that attempt 
to leverage the beneflts of both of the above modes of 
operation^^''^'''^^. On the other hand, those interested 
in implementing AQO strictly use the second mode of 
operation cited above. This, very naturally, leads to 
some interesting properties: First and foremost, in the 
coherent regime at e = 0, the groundstate maps onto 
Is) = (|0) + |1)) /V2, which implies that it is a super- 
position state with a fixed phase between components 
in the logical basis. Second, the logical basis is not co- 
incident with the energy eigenbasis, except in the ex- 
treme limit e/Aq ^ 1. As such, the qubit should not be 
viewed as an otherwise free spin-1/2 in a magnetic field, 
rather it maps onto an Ising spin subjected to a magnetic 
field with both a longitudinal (B^ — s- e) and a transverse 
{Bx Aq) component^^. In this case, it is the competi- 



tion between e and A^ which dictates the relative ampli- 
tudes of II) and It) in the groundstate wavefunction \g), 
thereby enabling logical operations that make no explicit 
use of the excited state |e). This latter mode of oper- 
ation of the flux qubit has connections to the fields of 
quantum magnetism^ and optimization theory Inter- 
estingly, systems of coupled flux qubits that are operated 
in this mode bear considerable resemblance to Feynman's 
original vision of how to build a quantum computer^^. 

While much seminal work has been done on single junc- 
tion and the related 3-Josephson junction rf-SQUID flux 

qubit20.21.22.23,24.25.26.27,28.29.30.31^ it haS been rCCOguized 

that such devices would be impractical in a large scale 
quantum information processor as their properties are 
exceptionally sensitive to fabrication variations. In par- 
ticular, in the regime Eji 3> hujp, Aq oc exp^—hiUp/ Eji). 
Thus, it would be unrealistic to expect a large scale pro- 
cessor involving a multitude of such devices to yield from 
even the best superconducting fabrication facility. More- 
over, implementation of AQO requires the ability to ac- 
tively tune Aq from being the dominant energy scale in 
the qubit to being essentially negligible during the course 
of a computation. Thus the single junction rf-SQUID flux 
qubit is of limited practical utility and has passed out of 
favor as a prototype qubit. 

The next step in the evolution of the single junction 
flux qubit and related variants was the compound Joseph- 
son junction (CJJ) rf-SQUID, as depicted in Fig. lb. 
This device was first reported upon by Han, Lapointe 
and Lukens'^^ and was the first type of fiux qubit to dis- 
play signatures of quantum superposition of macroscopic 
states'^'^. The CJJ rf-SQUID has been used by other re- 
search groups"'^'^"'^'*'^^ and a related 4-Josephson junction 
device has been proposed^o-^i. The CJJ rf-SQUID flux 
qubit and related variants have reappeared in a gradio- 
metric configuration in more recent history^^''^^''^^. Here, 
the single junction of Fig. la has been replaced by a flux 
biased dc-SQUID of inductance Lcjj that allows one to 
tune the critical current of the rf-SQUID in-situ. Let the 
applied flux threading this structure be denoted by <i>^jj . 
It is shown in Appendix A that the Hamiltonian for this 
system can be written as 
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where the sum is over n g {9,cjj}, Cq = C1+C2 
I/C1-I-I/C2 and Lq = Lbody + icjj/4. The 2-dimensional 
potential energy in Hamiltonian (4a) is characterized by 
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= 27r^ = - arctan ( ^ tan ((/3cjj/2) ) ; (4c) 
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Note that if cos (95cjj/2) < 0, then /3crr < in Haniiltonian 
(4a). This feature provides a natural means of shifting 
the qubit degeneracy point from (p^ = n, as in the single 
junction rf-SQUID case, to Pg « 0. It has been assumed 
in all that follows that this 7r-shifted mode of operation 
of the CCJ rf-SQUID has been invoked. 

Hamiltonian (4a) is similar to that of a single junc- 
tion rf-SQUID modulo the presence of a i^cjj-dependent 
tunnel barrier through PeB and an effective critical cur- 
l2- For Lcjj/Lq <C 1 it is reasonable to 



rent Ig = Ii 
assume that 



(Pcjj ~ ^7r<P;^jj/$o. Consequently, the CJJ 
rf-SQUID facilitates in-situ tuning of the tunneling en- 
ergy through While this is clearly desirable, one 
does pay for the additional flexibility by adding more 
complexity to the rf-SQUID design and thus more po- 
tential room for fabrication variations. The minimum 
achievable barrier height is ultimately limited by any so 
called junction asymmetry which leads to flnite In 
practice, for = {h - l2)/{h + h) < 0.05, this 
eff^ect is of little concern. However, a more insidious ef- 
fect of junction asymmetry can be seen via the change 
of variables pq — (p^ pq in Eq. (4a), namely an ap- 
parent $^jj-dependent flux offset: ff;^ - <i>°($^jj). 
If the purpose of the CJJ is to simply allow the experi- 
mentalist to target a particular A^, then the presence of 
<I>°(<I>|^jj) can be readily compensated via the application 
of a static flux offset. On the other hand, any mode of 
operation that explicitly requires altering Ag during the 
course of a quantum computation 
also require active compensation for 
a nonlinear crosstalk from to <&■ 
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.(.jj ^g. While it may be 
possible to approximate this effect as a linear crosstalk 
over a small range of if the junction asymmetry is 
small, one would nonetheless need to implement precise 
time-dependent flux bias compensation to utilize the CJJ 
rf-SQUID as a flux qubit in any quantum computation 
scheme. While this may be feasible in laboratory scale 
systems, it is by no means desirable nor practical on a 
large scale quantum information processor. 

A second problem with the CJJ rf-SQUID flux qubit 
is that one cannot homogenize the qubit parameters 1 | 
and Ag between a multitude of such devices that pos- 
sess different /3± over a broad range of ^cjj- While one 
can accomplish this task to a limited degree in a per- 
turbative manner about carefully chosen CJJ biases for 



each qubit ^, the equivalence of 



IP\ 



and Ao between 



those qubits will be approximate at best. Therefore, the 
CJJ rf-SQUID does not provide a convenient means of 
accommodating fabrication variations between multiple 
flux qubits in a large scale processor. 

Given that the CJJ rf-SQUID provides additional flex- 
ibility at a cost, it is by no means obvious that one can 
design a better rf-SQUID flux qubit by adding even more 
junctions. Specifically, it is desirable to have a device 
whose imperfections can be mitigated purely by the ap- 
plication of time-independent compensation signals. The 
novel rf-SQUID topology shown in Fig. Ic, hereafter re- 
ferred to as the compound-compound Josephson junction 



(CCJJ) rf-SQUID, satisfies this latter constraint. Here, 
each junction of the CJJ in Fig. lb has been replaced by a 
dc-SQUID, which will be referred to as left (L) and right 
(i?) minor loops, and will be subjected to external flux 
biases and respectively. The role of the CJJ loop 
in Fig. lb is now played by the CCJJ loop of inductance 
Lccjj which will be subjected to an external flux bias ^ccn ■ 
It is shown in Appendix B that if one chooses static val- 
ues of <I>^ and such that the net critical currents of 
the minor loops are equal, then it can be described by an 
effective two-dimensional Hamiltonian of the form 
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where the sum is over n G {5,ccjj}, Cg = Ci -\- C2 -\- 
C3 + C4, 1/Cccjj = l/(Ci + C2) + 1/(C3 + C4) and Lg = 
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energy in Hamiltonian (5a) is characterized by 
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where /3+($i,$^) = 27rig/g^($i, $^)/$o with 

I',{^1,^r) ^ (/1+/2)C0S (^!^)+(/3+4)cOS 



Given an appropriate choice of and the q and ccjj 
loops will possess apparent flux offsets of the form 
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where ^'^{r) is given by Eq. (B3c), which is purely a func- 
tion of ^L{R) junction critical currents. As such, 
the apparent flux offsets are independent of ^'ccjj- Un- 
der such conditions, we deem the CCJJ to be balanced. 
Given that the intended mode of operation is to hold 
and constant, then the offset phases p'j^ and p'j^ will 
also be constant. The result is that Hamiltonian (5a) 
for the CCJJ rf-SQUID becomes homologous to that of 
an ideal CJJ rf-SQUID [/3_ = in Eqs. (4b) and (4c)] 
with apparent static flux offsets. Such static offsets can 
readily be calibrated and compensated for in-situ using 
either analog control lines or on-chip programmable fltix 
sources^^. For typical device parameters and junction 
variability on the order of a few percent, these offsets will 
be ~ 1 ^ 10m$o- Equations 5a-5d with $° = $° -j = 
will be referred to hereafter as the ideal CCJJ rf-SQUID 
model. 

The second advantage of the CCJJ rf-SQUID is that 
one can readily accommodate for variations in criti- 
cal current between multiple flux qubits. Note that in 
Eq. (5b) that the maximum height of the tunnel barrier 



is governed by <I>fj) = + where 

is given by Eq. (B3c). One is free to choose any 
pair of (<i>|, <i>^) such that ^Li^f) = /3_r($^), as dictated 
by Eq. (B6). Consequently, (3+ = 2/3fl(4>^) in Eq. (5b). 
One can then choose <i>^, which then dictates $5^, so as 
to homogenize /3-|- between multiple flux qubits. The re- 
sults is a set of nominally uniform flux qubits where the 
particular choice of (<i>|^,<i>^) for each qubit merely re- 
sults in unique static flux offsets <i>° and ^Icjj ^^r each 
device. 

To summarize up to this point, the CCJJ rf-SQUID is 
robust against Josephson junction fabrication variations 
both within an individual rf-SQUID and between a plu- 
rality of such devices. The variations can be effectively 
tuned out purely through the application of static flux bi- 
ases, which is of considerable advantage when envisioning 
the implementation of large scale quantum information 
processors that use flux qubits. 



B. L-tuner 

The purpose of the CCJJ structure was to provide a 
means of coming to terms with fabrication variations in 
Josephson junctions both within individual flux qubits 
and between sets of such devices. However, junctions are 
not the only key parameter that may vary between de- 
vices, nor are fabrication variations responsible for all of 
the potential variation. In particular, it has been exper- 
imentally demonstrated that the inductance of a qubit 
Lq that is connected to other qubits via tunable mutual 
inductances is a function of the coupler configuration'*^. 
Let the bare inductance of the qubit in the presence of 
no couplers be represented by and the mutual induc- 
tance between the qubit and coupler i be represented by 
Mco,i- If the coupler possesses a first order susceptibihty 
Xi, as defined in Ref. 42, then the net inductance of the 
qubit can be expressed as 



(6) 



Given that qubit properties such as can be exponen- 
tially sensitive to variations in Lq, then it is undesirable 
to have variations in Lq between multiple flux qubits or 
to have Lq change during operation. This could have 
a deleterious impact upon AQO in which it is typically 
assumed that all qubits are identical and they are in- 
tended to be annealed in unison^. From the perspective 
of GMQC, one could very well attempt to compensate for 
such effects in a CJJ or CCJJ rf-SQUID flux qubit by ad- 
justing the tunnel barrier height to hold constant, but 
doing so alters | | , which then alters the coupling of the 
qubit to radiative sources, thus demanding further com- 
pensation. Consequently, it also makes sense from the 
perspective of GMQC that one find a means of render- 
ing Lq uniform between multiple qubits and insensitive 
to the settings of inductive coupling elements. 
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FIG. 3: (color online) A CCJJ rf-SQUID with L-tuner con- 
nected to multiple tunable inductive couplers via transformers 
with mutual inductances Mco,i and possessing susceptibilities 
Xi- The L-tuner is controlled via the external flux bias "1>i,t 



In order to compensate for variations in Lq, we have 
inserted a tunable Josephson inductance® into the CCJJ 
rf-SQUID body, as depicted in Fig. 3. We refer to this 
element as an inductance (L)-tuner. This relatively sim- 
ple element comprises a dc-SQUID whose critical current 
vastly exceeds that of the CCJJ structure, thus ensuring 
negligible phase drop across the L-tuner. Assuming that 
the inductance of the L-tuner wiring is negligible, the 
L-tuner modifies Eq. (6) in the following manner: 
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where Ljq = <i>o/27r/£2^, /£y is the net critical current of 
the two junctions in the L-tuner and *&f,T i^ externally 
applied flux bias threading the L-tuner loop. For modest 
flux biases such that L£y cos(7r<I>|j./<i>o) '> Iq, Eq. (7) is 
a reliable model of the physics of the L-tuner. 

Given that the L-tuner is only capable of augment- 
ing Lq, one can only choose to target Lq > Lq — 

J2i ^^co iXt^^ + Ljo, where xt^^ is the maximum anti- 
ferromagnetic (AFM) susceptibility of intcr-qubit coupler 
i. In practice, we choose to restrict operation of the cou- 
plers to the range —xt^^ < Xi < xt^^ such that the 
maximum qubit inductance that will be encountered is 
Lq> L° + Ei M^o,tXt^^ + Ljq. We then choose to pre- 
bias ^'-"^ each qubit to match the maximum realized 
amongst a set of flux qubits. Thereafter, one 
as couplers are adjusted by inverting 
Eq. (7) to solve for an appropriate value of Thus, 
the L-tuner provides a ready means of compensating for 
small variations in Lq between flux qubits and to hold 
Lq constant as inductive inter-qubit coupling elements 
are adjusted. 



III. DEVICE ARCHITECTURE, FABRICATION 
AND READOUT OPERATION 

To test the CCJJ rf-SQUID flux qubit, we fabricated 
a circuit containing 8 such devices with pairwise inter- 
actions mediated by a network of 16 in-situ tunable CJJ 
rf-SQUID inter-qubit couplers^^. Each qubit was also 
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FIG. 4: (color online) a) High level schematic of the analog 
devices on the device reported upon herein. Qubits are repre- 
sented as light grey elongated objects and denoted as go ... 97. 
One representative readout (RO), CCJJ and L-tuner (LT) 
each have been indicated in dashed boxes. Couplers (CO) 
are represented as dark objects located at the intersections of 
the qubit bodies, b) SEM of a cross-section of the fabrica- 
tion profile. Metal layers denoted as BASE, WIRA, WIRE 
and WIRC. Insulating layers labeled as Si02. Topmost in- 
sulator has not been planarized in this test structure, but is 
planarized in the full circuit process. An example via (VIA), 
Josephson junction (JUNC, AlO^^/Al) and resistor (RESI) are 
also indicated, c) Optical image of a portion of a device com- 
pleted up to WIRE. Portions of qubits qo . . .qs and the en- 
tirety of 54 are visible. 



coupled to its own dedicated quantum flux paranietron 
(QFP)-enabled readout''^. A high level schematic of the 
device architecture is shown in Fig. 4a. External flux bi- 
ases were provided to target devices using a sparse com- 
bination of analog current bias lines to facilitate device 
calibration and an array of single flux quantum (SFQ) 



based on-chip programmable control circuitry (PCC)''^. 

The device was fabricated from an oxidized Si wafer 
with Nb/Al/Al203/Nb trilayer junctions and four Nb 
wiring layers separated by planarized plasma enhanced 
chemical vapor deposited Si02. A scanning electron mi- 
crograph of the process cross-section is shown in Fig. 4b. 
The Nb metal layers have been labeled as BASE, WIRA, 
WIRB and WIRC. The flux qubit wiring was primarily 
located in WIRB and consisted of 2 /^m wide leads ar- 
ranged as an approximately 900 /im long differential mi- 
crostrip located 200 nm above a groundplane in WIRA. 
CJJ rf-SQUID coupler wiring was primarily located in 
WIRC, stacked on top of the qubit wiring to provide 
inductive coupling. PCC flux storage loops were im- 
plemented as stacked spirals of 13-20 turns of 0.25 /im 
wide wiring with 0.25 fj,m separation in BASE and WIRA 
(WIRB). Stored flux was picked up by one-turn washers 
in WIRB (WIRA) and fed into transformers for flux- 
biasing devices. External control lines were mostly lo- 
cated in BASE and WIRA. All of these control elements 
resided below a groundplane in WIRC. The groundplane 
under the qubits and over the PCC/external control lines 
were electrically connected using extended vias in WIRB 
so as to form a nearly continuous superconducting shield 
between the analog devices on top and the bias circuitry 
below. To provide biases to target devices with mini- 
mal parasitic crosstalk, transformers for biasing qubits, 
couplers, QFPs and dc-SQUIDs using bias lines and/or 
PCC elements were enclosed in superconducting boxes 
with BASE and WIRC forming the top and bottom, re- 
spectively, and vertical walls formed by extended vias in 
WIRA and WIRB. Minimal sized openings were placed 
in the vertical walls through which the bias and target 
device wiring passed at opposing ends of each box. 

An optical image of a portion of a device completed 
up to WIRB is shown in Fig. 4c. Qubits are visible as 
elongated objects, WIRB PCC spirals are visible as dark 
rectangles and WIRB washers are visible as light rect- 
angles with slits. Note that the extrema of the CCJJ 
rf-SQUID qubits are terminated in unused transformers. 
These latter elements allow this 8-qubit unit cell to be 
tiled in a larger device with additional inter-qubit CJJ 
couplers providing the connections between unit cells. 

We have studied the properties of all 8 CCJJ rf-SQUID 
flux qubits on this chip in detail and report upon one such 
device herein. To clearly establish the lingua franca of 
our work, we have depicted a portion of the multi-qubit 
circuit in Fig. 5a. Canonical representations of the exter- 
nal flux biases needed to operate a qubit, a coupler and 
a QFP-enabled readout are labeled on the diagram. The 
fluxes ^fj, <i>f^y and were only ever subjected to 
dc levels in our experiments that were controlled by the 
PCC. The remaining fluxes and readout current biases 
were driven by a custom-built 128 channel room temper- 
ature current source. The mutual inductance between 
qubit and QFP (Af,_qfp), between QFP and dc-SQUID 
(Afqfp.i-o), qubit and coupler {Mco,i) and <i>^Q-dependent 
inter-qubit mutual inductance (Mcff) have also been in- 
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FIG. 5; (color online) a) Schematic representation of a por- 
tion of the circuit reported upon herein. Canonical repre- 
sentations of all externally controlled flux biases "1?^, readout 
current bias iro and key mutual inductances Ma are indi- 
cated, b) Depiction of latching readout waveform sequence, 
c) Example QFP state population measurement as a function 
of the dc level $qfp with no qubit signal. Data have been fit 
to Eq. (8). 



dicated. Further details concerning cryogenics, magnetic 
shielding and signal filtering have been discussed in pre- 
vious publications'*^'^^'"'^''*'*. 

Since much of what follows depends upon a clear un- 
derstanding of our novel QFP-enabled readout mecha- 
nism, we present a brief review of its operation herein. 
The flux and readout current waveform sequence in- 
volved in a single-shot readout is depicted in Fig. 5b. 
Much like the CJJ qubit^'*, the QFP can be adiabat- 
ically annealed from a state with a monostable poten- 
tial (<&fatch — — *I'o/2) to a state with a bistable poten- 
tial (<&fatch = ^'^o) that supports two counter-circulating 
persistent current states. The matter of which persistent 
current state prevails at the end of an annealing cycle 
depends upon the sum of ^^^^ and any signal from the 
qubit mediated via A/q_qfp. The state of the QFP is then 
determined with high fidelity using a synchronized flux 
pulse and current bias ramp applied to the dc-SQUID. 
The readout process was typically completed within a 
repetition time ircp < 50 /iS. 

An example trace of the population of one of the QFP 
persistent current states Pqfp versus *&qfp, obtained using 
the latching sequence depicted in Fig. 5b, is shown in 
Fig. 5c. This trace was obtained with the qubit potential 
held monostable (*i>ccjj — ~^o/2) such that it presented 



minimal flux to the QFP and would therefore not influ- 
ence -Pqfp. The data have been flt to the phenomenolog- 
ical form 



1 — tanh 



qfp 



qfp 



(8) 



with width w ~ 0.18m<i>o for the trace shown therein. 



When biased with constant $ 



qfp 



$0 

qfp' 



which we refer 



to as the QFP degeneracy point, this transition in the 
population statistics can be used as a highly nonlinear 
flux amplifier for sensing the state of the qubit. Given 
that Mg_qfp = 6.28 ± 0.01 pH for the devices reported 
upon herein and that typical qubit persistent currents 
in the presence of negligible tunneling > 1 iiA, then 
the net fiux presented by a qubit was 2Afq^qfp \ Ip\ > 
6m<i>o, which far exceeded w. By this means one can 
achieve the very high qubit state readout fidelity reported 
in Ref. 43. On the other hand, the QFP can be used as a 
linearized flux sensor by engaging <i>qfp in a feedback loop 
and actively tracking 'i'qfp. This latter mode of operation 
has been used extensively in obtaining many of the results 
presented herein. 



IV. CCJJ RF-SQUID CHARACTERIZATION 

The purpose of this section is to present measurements 
that characterize the CCJJ, P-tuner and capacitance of 
a CCJJ rf-SQUID. All measurements shown herein have 
been made with a set of standard bias conditions given 
by $1 = 98.4m$o, '^R = -89.3 m$o, *lt = 0.344 $o 
and all inter-qubit couplers tuned to provide Meff = 0, 
unless indicated otherwise. The logic behind this partic- 
ular choice of bias conditions will be explained in what 
follows. This section will begin with a description of 
the experimental methods for extracting Lg and from 
persistent current measurements. Thereafter, data that 
demonstrate the performance of the CCJJ and i-tuner 
will be presented. Finally, this section will conclude with 
the determination of Cq from macroscopic resonant tun- 
neling data. 



High Precision Persistent Current 
Measurements 



The most direct means of obtaining information re- 
garding a CCJJ rf-SQUID is to measure the persistent 



current | | 



as a function of <f>?^^;;. A reasonable flrst 



approach to measuring this quantity would be to sequen- 
tially prepare the qubit in one of its persistent current 
states and then the other, and use the QFP in feedback 
mode to measure the difference in flux sensed by the QFP, 
which equals 2Mg_qfp \lq\- A fundamental problem with 
this approach is that it is sensitive to low frequency (LF) 
flux noise^^ , which can alter the background flux experi- 
enced by the QFP between the sequential measurements. 
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FIG. 6: (color online) a) Low frequency flux noise rejecting 
qubit persistent current measurement sequence. Waveforms 
shown are appropriate for measuring j/^ ('J'ccjj)l for — $0 < 
^ccjj ^ 0. The $ccjj waveform can be offset by integer $0 
to measure the periodic behavior of this quantity. Typical 
repetition time is trep ~ 1 ms. b) Depiction of QFP transition 
and correlated changes in QFP population statistics for the 
two different qubit initializations. 



flux bias is engaged in a differential feedback mode in 
which it is pulsed in alternating directions by an amount 
5^„i about a mean level The two single-shot mea- 
surements yield binary results for the QFP state and 
the difference between the two binary results is recorded. 
Gathering a statistically large number of such differential 
measurements then yields a differential population mea- 
surement (5Pqfp. Conceptually, the measurement works 
in the manner depicted in Fig. 6b: the two different ini- 
tializations of the qubit move the QFP degeneracy point 
to some unknown levels <i>^ ± (5<i>5^, where <I>^ represents 
the true mean of the degeneracy points at any given in- 
stant in time and 2(5<I>5^ is the true difference in degener- 
acy points that is independent of time. Focusing on flux 
biases that are close to the degeneracy point, one can 
linearize Eq. (8): 



qfp,± 



1 

2w 



(9) 



Assuming that the rms LF flux noise <i>„ <C w and that 
one has reasonable initial guesses for (f>°j ± (5<i>^, then 
the use of the linear approximation should be justified. 
Applying ^'^^^ = ± and sufficient repetitions of 
the waveform pattern shown in Fig. 6a, the differential 
population will then be of the form 



5P, 



qfp 



qfp.+ 



qfpr 



(10) 



which is independent of and <i>^. Note that the above 
expression contains only two independent variables, w 
and and that SPqip is purely a linear function of 

S^rn- By samphng at three values of S^m, as depicted by 
the pairs of numbered points in Fig. 6b, the independent 
variables in Eq. (10) will be overconstrained, thus readily 
yielding (5<I>^„. One can then infer the qubit persistent 
current as follows: 



, , 2(5$° 

jp ^^^m 

' 2Afq_qfp 



M, 



9 -qfp 



(11) 



For a typical measurement with our apparatus, the act of 
locating a single QFP degeneracy point to within 20 /i<i>o 
takes on the order of 1 s, which means that two sequential 
measurements would only be immune to ffux noise below 
0.5 Hz. We have devised a LF ffux noise rejection scheme 
that takes advantage of the fact that such noise will gen- 
erate a correlated shift in the apparent degeneracy points 
if the sequential preparations of the qubit can be inter- 
leaved with single-shot measurements that are performed 
in rapid succession. If these measurements are performed 



with repetition time t,-, 



1 ms, then the measurements 



will be immune to ffux noise below ^ 1 kHz. 

A depiction of the LF ffux noise rejecting persistent 
current measurement sequence is shown in Fig. 6a. The 
waveforms comprise two concatenated blocks of sequen- 
tial annealing of the qubit to a target ^ccH presence 
of an alternating polarizing ffux bias ±<i>^ followed by 
latching and single-shot readout of the QFP. The QFP 



Example measurements of i^'j:^;;) are shown in 
Fig. 7. These data, for which 1.5 < \Pcs\ S 2.5, have 
been fft to the ideal CCJJ rf-SQUID model by ffnding 
the value of ipq = </5™™ for which the potential in Eq. (5a) 
is minimized: 



\P'\ 



27r 



(12) 



The best fft shown in Fig. 7 was obtained with ~ 
265.4 ± 1.0 pH, Lccjj = 26 ± IpH and = 3.103 ± 
0.003 ^A. For comparison, we had estimated Lq = 273 pH 
at the standard bias condition for ^^rp and Lccjj = 20 pH 
from design. 

In practice, we have found that the LF ffux noise reject- 
ing method of measuring 1 | effectively eliminates any 
observable 1// component in that measurement's noise 
power spectral density, to within statistical error. Fi- 
nally, it should be noted that the LF ffux noise rejecting 
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FIG. 7: (color online) Example measurements of |/q (<i'ccjj)| 



method is applicable to any measurement of a difference 
in flux sensed by a linearized detector. In what follows 
herein, we have made liberal use of this technique to cali- 
brate a variety of quantities in-situ using both QFPs and 
other qubits as flux detectors. 



B. CCJJ 
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FIG. 8: (color online) Schematic of low frequency noise re- 
jecting differential qubit degeneracy point measurement se- 
quence. The qubit is annealed with a $ccjj signal of opposing 
polarity in the two frames and the qubit flux bias is controlled 
via feedback. 



In this subsection, the CCJJ has been characterized as 
a function of <I>2 ^"^d <I>|j, with all other static flux biases 
set to the standard bias condition cited above. Referring 
to Eq. (B4c), it can be seen that the qubit degeneracy 
point $° is a function of '^ccii through 70 if the CCJJ 
has not been balanced. To accentuate this functional 
dependence, one can anneal the CCJJ rf-SQUID with 
<i>^j.jj waveforms of opposing polarity about a minimum 
in |/3off|, as found at ^^cy] — ~*J'o/2- The expectation is 
that the apparent qubit degeneracy points will be anti- 
symmetric about the mean given by setting 70 = in 
Eq. (B4c). The waveform sequence for performing a dif- 
ferential qubit degeneracy point measurement is depicted 
in Fig. 8. In this case, the QFP is used as a latching 
readout and the qubit acts as the linearized detector of 
its own apparent annealing polarization-dependent flux 
offset. As with the |/^| measurement described above, 
this LF flux noise rejecting procedure returns a differ- 
ence in apparent flux sensed by the qubit and not the 
absolute flux offsets. 

To find balanced pairs of ((f>2i^fl:) i^i practice, we set 
<i>^ to a constant and used the LF flux noise rejecting 
procedure inside a software feedback loop that controlled 
<i>2 to null the difference in apparent degeneracy point to 
a precision of 20/Lt$o- Balanced pairs of (<&f,,<&|;) are 
plotted in Fig. 9a. These data have been fit to B6 using 
as a free parameter. The best fit shown in Fig. 9a 
was obtained with 1 - I3r^+/Pl,+ = (4.1 ± 0.3) x lO'^, 
which then indicates an approximately 0.4% asymmetry 



between the pairs of junctions in the L and R loops. 

A demonstration of how the CCJJ facilitates tuning of 
Iq is shown in Fig. 9b. Here, the measurable consequence 
of altering that was recorded was a change in |/|'| at 
([)■? .. = — $0- These data have been fit to the ideal CCJJ 
rf-SQUID model with the substitution 

p,{^i,<^i) = ncoB^-^^ (13) 

and using the values of iccjj and Lq obtained from fitting 
the data in Fig. 7, but treating /° as a free parameter. 
Here, <I>2 on the left side of Eq. (13) is a function of <i>^ 
per the CCJJ balancing condition Eq. (B6). The best fit 
was obtained with = 3.25 ±0.01 /j,A. This latter quan- 
tity agrees well with the designed critical current of four 
0.6 /im diameter junctions in parallel of 3.56 /iA. Thus, 
it is possible to target a desired 1^ by using Eq. (13) to 
select <i>^ and then Eq. (B6) to select <i>|_. The stan- 
dard bias conditions for <i>|^ and <i>i^ quoted previously 
were chosen so as to homogenize 1^ amongst the 8 CCJJ 
rf-SQUIDs on this particular chip. 



C. L- Tuner 

To characterize the L-tuner, we once again turned to 
measurements of |/^(<i>^^jj = — <i>o)|, but this time as a 
function of ^f^. Persistent current results were then 
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FIG. 9: (color online) a) Minor lobe balancing data and fit 
to Eq. (B6). The standard bias conditions for $f, and 3>fj are 
indicated by dashed lines, b) = — 'I'o)! versus 3>^, 

where has been chosen using Eq. (B6). The data have 
been fit to the ideal CCJJ rf-SQUID model. The standard 
bias condition for $fj and the resultant — — $o)| are 

indicated by dashed lines. 



used to infer SLq = L^(<i>|y) 



qi^LT — 0) using the 
and held con- 



ideal CCJJ rf-SQUID model with Lccjj 
stant and treating Lq as a free parameter. The experi 
mental results are plotted in Fig. 10a and have been fit 
to 



6L„ 



L 



JO 



COS (7r$f,y/$o) 



and the best fit was obtained with Ljq = 
Modeling this latter parameter as Lgo 



(14) 

19.60±0.04pH. 

= <i>o/27r/£j,, we 



estimate I^rp = 16.79 ± 0.04 /lA, which is close to the de- 
sign value of 16.94 /iA. The standard bias condition for 
<i>LT was chosen so as to homogenize Lq amongst the 8 
CCJJ rf-SQUID flux qubits on this chip and to provide 
adequate bipolar range to accommodate inter-qubit cou- 
pler operation. 

To demonstrate the use of the L-tuner, we have probed 
a worst-case scenario in which four CJJ rf-SQUID cou- 



FIG. 10: (color online) a) L-tuner calibration and fit to 
Eq. (14). The standard bias condition for $£x a-nd the resul- 
tant SLq are indicated by dashed lines, b) Observed change 
in maximum qubit persistent current with and without active 
L-tuner compensation and predictions for both cases. 



piers connected to the CCJJ rf-SQUID in question are 
tuned in unison. Each of the couplers had been in- 
dependently calibrated per the procedures described in 
Ref. 42, from which we obtained Mco,i « 15.8 pH and 
Xii^co) (* S {1,2,3,4}). Each of these devices pro- 
vided a maximum AFM inter-qubit mutual inductance 
Mafm = Afco,iXAFM ~ 1.56 pH, from which one can es- 
timate XAFM ~ 6.3 nH~^. Measurements of |J^| with 
and without active L-tuner compensation as a function 
of coupler bias <I>^q, as applied to all four couplers si- 
multaneously, are presented in Fig. 10b. The predictions 
from the ideal CCJJ rf-SQUID model, obtained by us- 
ing Lq = 265.4 pH (with compensation) and Lq obtained 
from Eq. (6) (without compensation), are also shown. 
Note that the two data sets and predictions all agree 
to within experimental error at <I>^q — 0.5<i>o, which 
corresponds to the all zero coupling state {Mcs ~ 0). 
The experimental results obtained without L-tuner com- 
pensation agree reasonably well with the predicted <i>co" 
dependence. As compared to the case without compensa- 
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tion, it can be seen that the measured \ show consider- 
ably less <i>|fQ-dependence when L-tuner compensation is 
provided. However, the data suggest a small systematic 
deviation from the inductance models Eqs. (6) and (7). 
At <i>^^jj — — for which it is estimated that PcS ~ 2.43, 
|/^| oc l/Lq. Given that the data for the case without 
compensation are below the model, then it appears that 
we have slightly underestimated the change in Lq . Conse- 
quently, we have provided insufficient ballast inductance 
when the L-tuner compensation was activated. 



D. rf-SQUID Capacitance 

Since 1^ and Lq directly impact the CCJJ rf-SQUID 
potential in Hamiltonian (5a), it was possible to in- 
fer CCJJ and L-tuner properties from measurements of 
the groundstate persistent current. In contrast, the rf- 
SQUID capacitance Cq appears in the kinetic term in 
Hamiltonian (5a). Consequently, one must turn to al- 
ternate experimental methods that invoke excitations of 
the CCJJ rf-SQUID in order to characterize Cq. One 
such method is to probe macroscopic resonant tunnel- 
ing (MRT) between the lowest lying state in one well 
into either the lowest order [LO, n = 0] state or into 
a higher order [HO, n > 0] state in the opposing well 
of the rf-SQUID double well potential^^. The spac- 
ing of successive HOMRT peaks as a function of rf- 
SQUID flux bias ^q will be particularly sensitive to 
Cq. HOMRT has been observed in many different rf- 
SQUIDs and is a well established quantum mechanical 
phenomenon^^^'*^''*^. LOMRT proved to be more difficult 
to observe in practice and was only reported upon rel- 
atively recently in the literature'*'*. We refer the reader 
to this latter reference for the experimental method for 
measuring MRT rates. 

Measurements of the initial decay rate T = dPi / dt\t=o 
versus <i>^ are shown in Fig. 11a with the order of the 
target level n as indicated. The maximum observable T 
was imposed by the bandwidth of the apparatus, which 
was ^ 5 MHz. The minimum observable F was dictated 
by experimental run time constraints. In order to ob- 
serve many HO resonant peaks within our experimental 
bandwidth we have successively raised the tunnel barrier 
height in roughly equal intervals by tuning the target 
^ccjj ■ The result is a cascade of resonant peaks atop a 
monotonic background. 

The authors of Ref. 46 attempted to fit their HOMRT 
data to a sum of gaussian broadened lorentzian peaks. It 
was found that they could obtain satisfactory fits within 
the vicinity of the tops of the resonant features but that 
the model was unable to correctly describe the valleys 
between peaks. We had reached the same conclusion with 
the very same model as applied to our data. However, it 
was empirically observed that we could obtain excellent 
fits to all of the data by using a model composed of a sum 
of purely gaussian peaks plus a background that varies 
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FIG. 11: (color online) a) HOMRT peaks fitted to Eq. (15). 
Data shown are for "l>?cjj/*o = -0.6677, -0.6735, -0.6793, 
-0.6851, -0.6911 and -0.6970, from left to right, respec- 
tively. Number of levels in target well n as indicated, b) Best 
fit Gaussian width parameter W„ as a function of n. c) Best 
fit peak position as a function of n. 



exponentially with 

where e = 2 |/|'| <i>^. These fits are shown in Fig. 11a. A 
summary of the gaussian width parameter Wn in Fig. lib 
is shown solely for informational purposes. We will re- 
frain from speculating why there is no trace of lorentzian 
lineshapes or on the origins of the exponential back- 
ground herein, but rather defer a detailed examination 
of HOMRT to a future publication. 

For the purposes of this article, the key results to 
take from the fits shown in Fig. 11a are the positions 
of the resonant peaks, as plotted in Fig. 11c. These 
results indicate that the peak spacing is very uniform: 
'5<i>MRT = 1.55±0.01 m<i>o. One can compare <5<i>MRT with 
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the predictions of the ideal CCJJ rf-SQUID model using 
the previously calibrated Lg = 265.4 pH, Lccjj = 26 pH 
and = 3.103 fiA with Cq treated as a free parameter. 
From such a comparison, we estimate Cq = 190 ± 2fF. 

The relatively large value of Cq quoted above can be 
reconciled with the CCJJ rf-SQUID design by noting 
that, unlike other rf-SQUID flux qubits reported upon in 
the literature, our qubit body resides proximal to a su- 
perconducting groundplane so as to minimize crosstalk. 
In this case, the qubit wiring can be viewed as a differ- 
ential transmission line of length £/2 ~ 900 /um, where 
£ is the total length of qubit wiring, with the effec- 
tive Josephson junction and a short on opposing ends. 
The transmission line will present an impedance of the 
form Z{uj) = —jZo ta.nh{uj£/2i') to the effective Joseph- 
son junction, with the phase velocity f = I/^/LqCq 
defined by the differential inductance per unit length 
Lq ^ 0.26pH//xm and capacitance per unit length Co ~ 
0.18 fF /^m, as estimated from design. If the separation 
between differential leads is greater than the distance to 
the groundplane, then £/2i/ « ■\/.'jbodyC'body/4, where 
Cbody ~ 640 fF is the total capacitance of the qubit 
wiring to ground. Thus, one can model the high fre- 
quency behavior of the shorted differential transmission 
line as an inductance Lbody and a capacitance Cbody /4 
connected in parallel with the CCJJ. Taking a reason- 
able estimated value of 40 fF/ /im^ for the capacitance 
per unit area of a Josephson junction, one can esti- 
mate the total capacitance of four 0.6 fxm diameter junc- 
tions in parallel to be Cj ~ 45 fF. Thus we estimate 
Cq = Cj + Cbody /4 ^ 205 fF, which is in reasonable 
agreement with the best fit value of Cq quoted above. 

With all of the controls of the CCJJ rf-SQUID hav- 
ing been demonstrated, we reach the first key conclusion 
of this article: The CCJJ rf-SQUID is a robust device 
in that parametric variations, both within an individual 
device and between a multitude of such devices, can be 
accounted for using purely static flux biases. These bi- 
ases have been applied to all 8 CCJJ rf-SQUIDs on this 
particular chip using a truly scalable architecture involv- 
ing on-chip flux sources that are programmed by only a 
small number of address lines*^. 



V. QUBIT PROPERTIES 

The purpose of the CCJJ rf-SQUID is to provide an 
as ideal as possible flux qubit ^. By this statement, it is 
meant that the physics of the two lowest lying states of 
the device can be described by an effective Hamiltonian 
of the form Eq. (3) with e = 2 (<I>;^ - $^), being 
the magnitude of the persistent current that flows about 
the inductive loop when the device is biased hard to one 
side, being a static flux offset and Aq representing the 
tunneling energy between the lowest lying states when bi- 
ased at its degeneracy point = Thus, and 
are the deflning properties of a flux qubit, regardless of its 
topology^. Given the complexity of a six junction device 



with flve closed superconducting loops, it is quite justifl- 
able to question whether the CCJJ rf-SQUID constitutes 
a qubit. These concerns will be directly addressed herein 
by demonstrating that measured and Ag agree with 
the predictions of the quantum mechanical Hamiltonian 
(5a) given independently calibrated values of Lq, iccjj) 
and Cq. 

Before proceeding, it is worth providing some con- 
text in regards to the choice of experimental methods 
that have been described below. For those researchers 
attempting to implement GMQC using resonant elec- 
tromagnetic flelds to prepare states and mediate inter- 
actions between qubits, experiments that involve high 
frequency pulse sequences to drive excitations in the 
qubit (such as Rabi oscillations^^, Ramsey fringes^^''^° 
and spin-echo^^'^°'^^) are the natural modality for study- 
ing quantum effects. Such experiments are convenient in 
this case as the methods can be viewed as basic gate oper- 
ations within this intended mode of operation. However, 
such methods are not the exclusive means of characteriz- 
ing quantum resources. For those who wish to use precise 
dc pulses to implement GMQC or whose interests lie in 
developing hardware for AQO, it is far more convenient 
to have a set of tools for characterizing quantum mechan- 
ical properties that require only low bandwidth bias con- 
trols. Such methods, some appropriate in the coherent 
regime'*^'^^ and others in the incoherent regime^^''*^'^^, 
have been reported in the literature. We have made use 
of such low frequency methods as our apparatuses typ- 
ically possess 128 low bandwidth bias lines to facilitate 
the adiabatic manipulation of a large number of devices. 

One possible means of probing quantum mechanical 
tunneling between the two lowest lying states of a CCJJ 
rf-SQUID is via MRT^^. Example LOMRT decay rate 
data are shown in Fig. 12a. We show results for both 
initiahzations, \l) and |t), and fits to gaussian peaks, as 
detailed in Ref. 44: 

r($-) /^A2 

A summary of the fit parameters Cp and W versus ^ccjj 
is shown in Fig. 12b. We also provide estimates of the 
device temperature using the formula 

fefilMRT = • (17) 

As expected, Tmrt shows no discernible ^ccjj' 
dependence and is scattered about a mean value of 
53 ± 2 mK. A summary of A^ versus ^ccjj ^^^^ shown 
in conjunction with more experimental results at the end 
of this section. For further details concerning LOMRT, 
the reader is directed to Ref. 44. 

A second possible means of probing A^ is via a Landau- 
Zener experiment^". In principle, this method should be 
applicable in both the coherent and incoherent regime. In 
practice, we have found it only possible to probe the de- 
vice to modestly larger A^ than we can reach via LOMRT 
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FIG. 12: (color online) a) Example LOMRT peaks fitted to 
Eq. (16). Data shown are for $?cjj/*o = -0.6621, -0.6642 
and —0.6663, from top to bottom, respectively. Data from 
the qubit initialized in |J.) (|t)) are indicated by solid (hol- 
low) points, b) Energy scales obtained from fitting multiple 
LOMRT traces. 



purely due to the low bandwidth of our bias lines. Results 
from such experiments on the CCJJ rf-SQUID flux qubit 
will be summarized at the end of this section. We see no 
fundamental limitation that would prevent others with 
higher bandwidth apparatuses to explore the physics of 
the CJJ or CCJJ flux qubit at the crossover between the 
coherent and incoherent regimes using the Landau-Zener 
method. 

In order to probe the qubit tunnel splitting in the co- 
herent regime using low bandwidth bias lines, we have 
developed a new experimental procedure for sensing the 
expectation value of the qubit persistent current, simi- 
lar in spirit to other techniques already reported in the 
literature*^. An unfortunate consequence of the choice 
of design parameters for our high fidelity QFP-enabled 
readout scheme is that the QFP is relatively strongly cou- 
pled to the qubit, thus limiting its utility as a detector 
when the qubit tunnel barrier is suppressed. One can cir- 
cumvent this problem within our device architecture by 



tuning an inter-qubit coupler to a finite inductance and 
using a second qubit as a latching sensor, in much the 
same manner as a QFP. Consider two flux qubits coupled 
via a mutual inductance M^s- The system Hamiltonian 
can then be modeled as 



ie{qM} 



f^+A,a^A+Ja<f^af^ , (15 



where J = Meg\I^\\I^\. Let qubit q be the flux source 
and qubit d serve the role of the detector whose tun- 
nel barrier is adiabatically raised during the course of a 
measurement, just as in a QFP single shot measurement 
depicted in Fig. 5. In the limit A^; — )■ one can write 
analytic expressions for the dispersion of the four lowest 
energies of Hamiltonian (18): 



Ei± 

E2± 



2<^d 



(19) 



As with the QFP, let the flux bias of the detector qubit be 
engaged in a feedback loop to track its degeneracy point 
where P^^i = 1/2. Assuming Boltzmann statistics for the 
thermal occupation of the four levels given by Eq. (19), 
this condition is met when 



,-E2-/kBT _|_ ^-E2+/kBT 



y 



e{i±,2±} 



g-E^/kBT 



(20) 



Setting Pd^i = 1/2 in Eq. (20) and solving for £2 then 
yields an analytic formula for the balancing condition: 



fid = 



F{+)-F{-) 



kgTlll 



1 + p-F{ + )/kBT 
1 + e-F(-)/feBT 



(21) 



Fi±)^d{e,±2jr + Al. 



While Eq. (21) may look unfamiliar, it readily reduces 
to an intuitive result in the limit of small coupling J <C 
Ai and T ^ 0: 



Mcff (5I Iff) 



(22) 



where \g) denotes the groundstate of the source qubit 



and IP = (7 



(9) 



is the source qubit persistent current 
operator. Thus Eq. (21) is an expression for the expecta- 
tion value of the source qubit 's groundstate persistent 
current in the presence of backaction from the detec- 
tor and flnite temperature. Setting = 2|/f|<i>f and 
rearranging then gives an expression for the flux bias 
of the detector qubit as a function of flux bias applied 
to the source qubit. Given independent calibrations of 
Afcff = 1-56 ± 0.01 pH for a particular coupler set to 
= on this chip, T = 54 ± 3 mK from LOMRT 
flts and = 1.25 ± 0.02 ^A at the CCJJ bias where 
the LOMRT rate approaches the bandwidth of our bias 
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FIG. 13: (color online) Example coupled flux trace taken at 
^ccjj = —0.6513 $0 used to extract large A in the coherent 
regime. 



lines, one can then envision tracing out (f>^ versus <i>^ and 
fitting to Eq. (21) to extract the source qubit parameters 
\IP\ and A, . 

An example <i>^ versus <i>^ data set for source CCJJ flux 
bias <i>|?j.jj = -0.6513 <i>o is shown in Fig. 13. The sohd 
curve in this plot corresponds to a fit to Eq. (21) with 
a small background slope that we denote as x- We have 
confirmed from the ideal CCJJ rf-SQUID model that x is 
due to the diamagnetic response of the source rf-SQUID 
to changing <I>^. This feature becomes more pronounced 
with increasing Cq and is peaked at the value of <&ccjj for 
which the source qubit potential becomes monostable, 
PeS = 1- Nonetheless, the model also indicates that x in 
no way modifies the dynamics of the rf-SQUID, thus the 
qubit model still applies. From fitting these particular 
data, we obtained = 0.72 ± 0.04 and Ag/h = 
2.64 ±0.24 GHz. 

In practice we have found it inefficient to take detailed 
traces of versus as this procedure is susceptible 
to corruption by LF flux noise in the detector qubit. As 
an alternative approach, we have adapted the LF flux 
noise rejecting procedures introduced in the last section 
of this article to measure a series of three differential 
flux levels in the detector qubit. The waveforms needed 
to accomplish this task are depicted in Fig. 14. Here, 
the dc-SQUID and QFP connected to the detector qubit 
are used in latching readout mode while the detector 
qubit is annealed in the presence of a differential flux bias 
± S^,n which is controlled via feedback. Meanwhile, 
the source qubit 's CCJJ bias is pulsed to an intermediate 
level — <i>o < ^ccjj < — *i'o/2 in the presence of an initial- 
ization flux bias ±<i>g. By choosing two appropriate pairs 
of levels ±$^, as indicated by the solid points 1± and 2± 
in Fig. 13, one can extract |/|'| and x from the two dif- 
ferential flux measurements. In order to extract A^, we 
then choose a pair of ±<i>^ in the centre of the trace, as 
indicated by the solid points 3±, from which we obtain 
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FIG. 14: (color online) Depiction of large measurement 
waveforms. The waveform sequence is similar to that of Fig. 6, 
albeit the source qubit's tunnel barrier is partially suppressed 
(— $o/2 < "I'ccjj < — 'I'o) and a second qubit (as opposed to a 
QFP) serves as the flux detector. 



the central slope d^^/d^^. Taking the first derivative of 
Eq. (21) and evaluating at <i>^ = yields 



d^ 
d<^^ 



X = 



: tanh 



{2JY 



(2J)^ + 
2kbT 



(23) 



Given independent estimates of all other parameters, one 
can then extract A^ from this final differential flux mea- 
surement . 

A summary of experimental values of the qubit pa- 
rameters 1 1^1 and Ag versus <i>^^jj is shown in Fig. 15. 
Here, we have taken A^ from LOMRT and Landau- 
Zener experiments in the incoherent regime and from 
the LF flux noise rejecting persistent current procedure 
discussed above in the coherent regime. The large gap 
between the three sets of measurements is due to two rea- 
sons: First, the relatively low bandwidth of our bias lines 
does not allow us to perform MRT or Landau-Zener mea- 
surements at higher Ag where the dynamics are faster. 
Second, while the coherent regime method worked for 
Ag > fcsT, it proved difficult to reliably extract A^ in 
the opposite limit. As such, we cannot make any pre- 
cise statements regarding the value of <i>^^jj which serves 
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15: (color online) a) Magnitude of the persistent cur- 

b) Tunneling energy A, be- 



\J-q I 



as a function of $ccjj • 



tween two lowest lying states of the CCJJ rf-SQUID as a 
function of $ccjjj characterized by macroscopic resonant 
tunneling [MRT] and Landau-Zener [LZ] in the incoherent 
regime and coupled groundstate persistent current {{g\ I'^ \g)) 
in the coherent regime. Solid curves are the predictions of the 
ideal CCJJ rf-SQUID model using independently calibrated 
Lq, Lccjj, Iq and Cq with no free parameters. 
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FIG. 16: (color online) Low frequency flux noise in the CCJJ 
rf-SQUID flux qubit. Data [points] have been fit to Eq. (24) 
[solid curve]. 



be identified as a flux qubit as the measured \lq\ and 
Ag agree with the predictions of a quantum mechanical 
Hamiltonian whose parameters were independently cali- 
brated. 



VI. NOISE 

With the identification of the CCJJ rf-SQUID as a 
flux qubit firmly established, we now turn to assessing 
the relative quality of this device in comparison to other 
flux qubits reported upon in the literature. In this sec- 
tion, we present measurements of the low frequency flux 
and critical current spectral noise densities, S^{f) and 
Si{f), respectively. Finally, we provide explicit links be- 
tween S^{f) and the free induction (Ramsey) decay time 
T2 that would be relevant were this flux qubit to be used 
as an element in a gate model quantum information pro- 
cessor. 



as the delineation between the coherent and incoherent 
regimes based upon the data shown in Fig. 15b. Reg- 
ulating the device at lower temperature would assist in 
extending the utility of the coherent regime method to 
lower Ag. On the other hand, given that Eq. (16) predicts 
that F cx Ag, one would have to augment the experimen- 
tal bandwidth by at least two orders of magnitude to 
gain one order of magnitude in Ag via either MRT or LZ 
experiments. 

The solid curves in Fig. 15 were generated with the 
ideal CCJJ rf-SQUID model using the independently cal- 
ibrated Lq = 265.4 pH, Lccjj = 26 pH, I^j = 3.103 ^iA and 
Cq = 190 fF. Note that there are no free parameters. It 
can be seen that the agreement between theory and ex- 
periment is quite reasonable. Thus we reach the second 
key conclusion of this article: The CCJJ rf-SQUID can 



A. Flux Noise 

Low frequency (1//) flux noise is ubiquitous in super- 
conducting devices and is considered a serious impedi- 
ment to the development of large scale solid state quan- 
tum information processors'*^. We have performed sys- 
tematic studies of this property using a large number of 
flux qubits of varying geometry^ ^ and, more recently, as a 
function of materials and fabrication parameters. These 
latter studies have aided in the reduction of the ampli- 
tude of 1// flux noise in our devices and will be the 
subject of a forthcoming publication. Using the methods 
described in Ref. 51, we have generated the one-sided 
flux noise power spectral density 5'$(/) shown in Fig. 16. 
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These data have been fit to the generic form 



Sif) 



(24) 



with best fit parameters a = 0.95 ± 0.05, ^/w^ = 9.7 ± 
0.5/i<I>o/\/Hz and amplitude A such that 5$(1 Hz) = 
l-3lo 5 /^'^'o/VHz. Thus wc reach the third key conclu- 
sion of this article: Wc have demonstrated that it is 
possible to achieve 1// flux noise levels with Nb wiring 
that are as low as the best Al wire qubits reported in 
the literature'^*''^"'^''*'''. Moreover, we have measured simi- 
lar spectra from a large number of identical flux qubits, 
both on the same and different chips, and can state with 
confidence that the 1// amplitude reported herein is re- 
producible. Given the experimentally observed geometric 
scaling of 5$(1 Hz) in Ref. 51 and the relatively large size 
of our flux qubit bodies, we consider the prospects of ob- 
serving even lower 1// noise in smaller flux qubits from 
our fabrication facility to be very promising. 



B. Critical Current Noise 

A second noise metric of note is the critical current 
noise spectral density Si{f). This quantity has been 
studied extensively and a detailed comparison of exper- 
imental results is presented in Ref. 52. A recent study 
of the temperature and geometric dependence of criti- 
cal current noise has been published in Ref. 53. Based 
upon Eq. (18) of Ref. 52, we estimate that the 1// crit- 
ical current noise from a single 0.6 /im diameter junc- 
tion, as found in the CCJJ rf -SQUID flux qubit, will 
have an amplitude such that -\/ 5/(1 Hz) <^ 0.2pA/\/Hz. 
Unfortunately, we were unable to directly measure criti- 
cal current noise in the flux qubit. While the QFP-enable 
readout provided high fldelity qubit state discrimination 
when qubits are fully annealed to ^^cn ~ ~*J'o, this read- 
out mechanism simply lacked the sensitivity required for 
performing high resolution critical current noise measure- 
ments. In lieu of a measurement of S'/(/) from a qubit, 
we have characterized this quantity for the dc-SQUID 
connected to the qubit in question. The dc-SQUID had 
two 0.6 jMQ. junctions connected in parallel. A time trace 
of the calibrated switching current /gw ~ Ic was obtained 
by repeating the waveform sequence depicted in Fig. 5b 
except with ^f^tch = — *I'o/2 at all time (QFP disabled, 
minimum persistent current) and = to provide 
minimum sensitivity to flux noise. Assuming that the 
critical current noise from each junction is uncorrelated, 
the best that we could establish was an upper bound 
of ■\/S'/(lHz) < 7pA/v^Hz for a single 0.6 /um diameter 
junction. 

Given the upper bound cited above for critical current 
noise from a single junction, we now turn to assessing the 
relative impact of this quantity upon the CCJJ rf-SQUID 
flux qubit. It is shown in Appendix B that fluctuations in 
the critical currents of the individual junctions of a CCJJ 



generate apparent flux noise in the flux qubit by modu- 
lating Inserting critical current fluctuations of mag- 
nitude 5Ic < 7pA/v^Iz and a mean junction critical cur- 
rent Ic = ~ 0.8 ^A into Eq. (BIO) yields qubit de- 
generacy point fluctuations ^ 0.1 /Lt$o/V^Hz- This 
final result is at least one order of magnitude smaller than 
the amplitude of 1// fiux noise inferred from the data in 
Fig. 16. As such, we consider the effects of critical current 
noise in the CCJJ rf-SQUID to be tolerable. 



C. Estimation of 

While measurements of noise power spectral densities 
are the most direct way of reporting upon and comparing 
between different qubits, our research group is frequently 
asked what is the dephasing time for our fiux qubits. 
The answer presumably depends very strongly upon bias 
settings, for recall that we have measured properties of 
the CCJJ rf-SQUID fiux qubit in both the coherent and 
incoherent regime. Given that our apparatuses contain 
only low bandwidth bias lines for enabling AQO, we are 
unable to measure dephasing within our own laboratory. 
Collaborative efforts to measure dephasing for our fiux 
qubits are in progress. In the meantime, we provide a 
rough estimate below for our fiux qubits if they were 
biased to the optimal point, <I>^ = <I>° based upon the 
measured S'$(/) and subjected to a free induction decay, 
or Ramsey fringe, experiment. Referring to Eq. (33a) of 
Ref. 54 and key results from Ref. 55, the mean squared 
phase noise for a flux qubit at the optimal point will be 
given by 



{€{t)) 



1 {m\) 

ffl 2A2 



^ ..^ sin^(^/t) .25i 



where S^2{ijj) represents the quadratic fiux noise spec- 
tral density and is the measurement cutoff frequency. 
Assuming that the first order spectral density S^{uj) = 
27rA^/a;, then S^2{w) can be written as 



^ /die-- /d.'^ /d."^ 

2tt J J CO' J LO" 



27r 

= 877^^^ 



In (w/wir) 



(26) 



where Wir = 27r/ir denotes an infrared cutoff of the 1// 

noise spectral density. Inserting Eq. (26) into Eq. (25) 
and rendering the integral dimensionless then yields: 

J mm*' 

(27) 

where /,nin = max [ /,„ /ir ] . We have numerically stud- 
ied the behavior of the integral in Eq. (27). In the very 
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long measurement time limit the integral is cut off by /ir 
and the integral varies as 1/t^, which then cancels the 
factor of in the numerator of Eq. (27). This means 
that the mean squared phase noise eventually reaches a 
finite limit. However, the more experimentally relevant 
limit is fm ^ /ir , for which we found empirically that 
the integral varies roughly as 5 x [In (/m/Zir)]^ over many 
orders of magnitude in the argument of the logarithm. In 
this latter limit the result is independent of t, so Eq. (27) 
can be rewritten as {4''^{t)) = /{T^)^, which then yields 
the following formula for : 



1 (2|/f|A)^ 
7rA2 



-1-1/2 



5 In (/„//;,) 



(28) 



Since fiux noise spectra seem to obey the 1// form 
down to at least 0.1 mHz and researchers are generally 
concerned with dephasing over times of order 1 /zs, then 
it is fair to consider fm/ fir ~ 10^*^. For a nominal value 
of such that the flux qubit is in the coherent regime, 
say — 0.652 $0) the qubit parameters are Ag/Zi « 2 GHz 
and w 0.7 /zA. Substituting these quantities into 
Eq. (27) then yields ^ 150 ns. This estimate of the 
dephasing time is comparable to that observed in consid- 
erably smaller flux qubits with comparable 1// flux noise 
levels30>3i. 



V'cjj = ^Pl- ^2 



(Alb) 



respectively. Furthermore, assume that the CJJ loop has 
an inductance Lqj that is divided symmetrically between 
the two paths. Using trigonometric relations, one can 
write a Hamiltonian for this system in terms of modes in 
the q and cjj loops that has the following form: 



H 



E 



2C, 



-Ugf3+ COS COS {ifg) 

+f/g/3_sin(^)sin(^g) ; 



27rL, (Ji ± h) 



(A2a) 



(A2b) 



where the sum is over n G {9, cjj}, Cq = C1+C2, l/Ccjj = 
I/C1+I/C2, Un = ($o/27r)VL„, Lq = Lbody+i-cjj/4and 
[i>o<Pn/27r, Qn] = ih. The Josephson potential energy 
of Hamiltonian (A2a) can be rearranged by defining an 
angle 9 such that tan^ = (/3_//3+) tan (iy9cjj/2). Further 
trigonometric manipulation then yields Eqs. (4a)-(4d). 



APPENDIX B: CCJJ RF-SQUID 



VII. CONCLUSIONS 

One can draw three key conclusions from the work pre- 
sented herein: First, the CCJJ rf-SQUID is a robust and 
scalable device in that it allows for in-situ correction for 
parametric variations in Josephson junction critical cur- 
rents and device inductance, both within and between 
fiux qubits using only static flux biases. Second, the mea- 
sured fiux qubit properties, namely the persistent current 
1 71* I and tunneling energy Aq, agree with the predictions 
of a quantum mechanical Hamiltonian whose parameters 
have been independently calibrated, thus justifying the 
identification of this device as a fiux qubit. Third, it 
has been experimentally demonstrated that the low fre- 
quency fiux noise in this all Nb wiring fiux qubit is com- 
parable to the best all Al wiring devices reported upon 
in the literature. Taken in summation, these three con- 
clusions represent a significant step forward in the devel- 
opment of useful large scale superconducting quantum 
information processors. 

We thank J. Hilton, P. Spear, A. Tcaciuc, F. Cioata, 
M. Amin, F. Brito, D. Averin, A. Kleinsasser and G. Ker- 
ber for useful discussions. Siyuan Han was supported in 
part by NSF Grant No. DMR-0325551. 



Following the same logic as for the CJJ rf-SQUID, one 
can define four orthogonal quantum mechanical degrees 
of freedom as follows: 

= ^1-^2; (Bla) 
(PR = ips-Lfi] (Bib) 

If I + 1^2 + 'PA 



V'ccjj = 'Pi — 'Pr = 
'Pq 



2 2 

Ve + Vr _ ^1 + ^2 + ^3 + 

2 ~ 4 



; (Blc) 
. (Bid) 



Using the same strategy as in Appendix A, one can use 
trigonometric identities to first express the Josephson po- 
tential in terms of the L and R loop modes: 

^ = Z.^ + Z.f^" 2 

n m 
-Uql3L+ COS (^) COS {iplj 

+Uql3L- sin ^'^^ sin {ipi) 

-Uq(3R+ cos cos {(fir) 

+Uqf3R_ sin sin {ip^) ; (B2a) 



APPENDIX A: CJJ RF-SQUID 

Let the qubit and cjj loop phases be defined as 

ifiq = {(fii + <p2) /2 , (Ala) 



_ 27rLg (7i(3) ±72(4)) 

fiL{R),± = 



(B2b) 



where the first sum is over n € {L, R, r} and the second 
sum is over closed inductive loops m G {L, R, ccjj, q}. As 
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before, each of the modes obey the commutation relation 
[^otPn/ZTT, Q„] = ih. Here, 1/Cl(r) = 1/Ci(3) + 1/C2(4), 
Q(r) = C'i(3) + C'2(4) and Um = ($o/27r)=^/i„. 

We have found it adequate for our work to assume 
that LL_ii/Lq <C 1; which then allows one to reduce 
the four dimensional system given in Hamiltonian (B2a) 
to two dimensions. Consequently, we will substitute 
fLiR) = Vl(r) ignore the L and R kinetic terms 
henceforth. Assuming that the inductance of the ccjj 
loop is divided equally between the two branches one can 
then write Lq = Lbody + ^ccjj / 4. With these approxima- 
tions and the 6 strategy presented in Appendix A, one 
can rearrange the Josephson potential terms to yield the 
following: 



n 



L{R) 



-UqiSL COS {ife - (fl) 
-UqPR COS {ipr - ^n) ; 

13 L{R),+ COS ( 1 



(B3a) 



(B3b) 




Pl(R)- ( Vl{R) 



/3_ 



l[r:),+ 



V 2 



<^L(fl) = -arctan(^^^|^tan(v^^(^)/2)) , (B3c) 

where the first sum is over n G {£, r} and the second 
sum is over m G {cc]],q\. The Josephson potential is 
given by a sum of two cosines, as encountered in the CJJ 
rf-SQUID derivation of Hamiltonian (A2a) from Hamil- 
tonian (1). These two terms can be rewritten in the same 
manner by defining 0± = /3i ± (Sr. The result, similar 
to Hamiltonian (A2a), can then be subjected to the 6 
strategy to yield 



n 



E 



01 

2C„ 



-f/,/3effC0S (</9q - (^°) , 



(B4a) 



where the sum is over n S ccjj} and the capacitances 
are defined as Cq = Ci + C2 + Cz + Ci and 1/C, 



+ C2) 
defined as 



CCJJ 



1/(C3 + Ci). The other parameters are 



/3eff = /3+cosh^ Wl + 



/3+ 



tan — 



(i) 



'fa = n ^ 70 ; 



7 = <^ccjj 



(B4b) 

(B4c) 
(B4d) 



7o = — arctan ( ^ tan(7/2) 

P± = Pl±Pr. 



(B4e) 



(B4f) 



Hamiltonian (B4a) inherits much of its complexity 
from junction asymmetry both within the minor loops, 
which gives rise to t/'L(_R)) and effective junction asymme- 
try between the minor loops, which gives rise to 70. For 
arbitrary external fiux biases and nominal spread injunc- 
tion critical current, the CCJJ rf-SQUID offers no obvi- 
ous advantage over the CJJ rf-SQUID. However, upon 
choosing biases and such that 



0L = 0R 



(B5) 



then /?_ = 0, and consequently 70 = 0. With these sub- 
stitutions, Hamiltonian (B4a) yields Hamiltonian (5a). 

Note that for 0l(r).- / Pl(r).+ ^ 1 and modest 

that the so-called CCJJ balancing condition given by 

Eqs. (B3b) and (B5) can be written approximately as 



13l,+ cos 



ft 
2 



fiR,+ COS 
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which, upon solving for ^p'f^ yields 
27r 



arccos 



^-'-cos^*« 



(3. 



$0 



(B6) 



It is possible for critical current noise to couple into 
the ipq degree of freedom in any compound junction rf- 
SQUID qubit via modulation of the junction asymmetry- 
dependent apparent qubit flux offset <I>°. In the case of 
the CCJJ rf-SQUID, all three quantities on the right side 
of Eq. (B4c) are ultimately related to the critical cur- 
rents of the individual junctions. Given typical junction 
parameter spreads from our fabrication facility. 



0(0.01) , 



Pl(r)- 




Il(3) 


- ^2(4) 


I3l{r),+ 




A (3) 


+ ^2(4) 



SO one can write an approximate expression for </?2,(k) 
using Eq. (B3c): 
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and for 70 using Eqs. (B3b), (B4d) and (B4e): 
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where Ic represents the mean critical current of a single 
junction. The CCJJ rf-SQUID is intended to be oper- 
ated with only small flux biases in the minor loops, thus 

cos {^^2^ ~ cos ~ 1- It is also reasonable to as- 

sume that 7 « feci) ^ corrections to tan(7/2) from 
V>\(ii) and from the effective two-dimensionality of the 
rf-SQUID potential will be very small. Inserting Eqs. B7 
and B8 into Eq. (B4c) then yields 



article, ^'^(^j/^'o ^ 0.1 and the device acts as a qubit for 
$^j.jj/$o ~ 0.65. For these flux biases, the magnitude of 
the terms within the square braces in Eq. (B9) are all of 
order 1. Therefore, for general flux bias conditions, the 
apparent qubit flux offset is roughly given by 
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"47r L 
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- tan 
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+ - tan 



tan 

2 



2 



— tan 



fk 
2 



tan 



2 

2 



,(B9) 



For the typical operating parameters described in this 



Assume that each junction experiences critical current 

fluctuations of magnitude Sic- If each junction's fluctua- 
tions axe independent, then the root mean square varia- 
tion of the qubit degeneracy point will be 



^5h 
27r I, 



(BIO) 



Thus, critical current fluctuations generate apparent flux 
noise in the CCJJ rf-SQUID flux qubit. 
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